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^ ■ 1 INTRODUCTION 
7-H ■ 

The role of globular clusters (GCs) in the chemical evolu- 
. ■ tion of galaxies will be better understood by studying car- 
bon stars (CSs) in the former. Only two types of CSs have 
?H \ abundances and kinematics typical for the Galactic halo : the 
^ giant and dwarf CH stars (|McClureHl985l ; lGreenlll99rj l. In 
the field most of these stars are binaries, where the primary 
owes its peculiar nature to mass transfer from its white- 
dwarf companion rather than to dredge-up from the interior 
jMcClure fc Woodsworthlll99ol l. There is no observational 
evidence for the exact nature of CH stars in GCs. The chem- 
ical differences between CH stars that are intrinsic asymp- 
totic giant bra nch (AGB) stars a nd those in binary systems 
are detailed bv lAbia etall (|2003l V 

Only three CH stars wer e found in GCs so far: two in 



ABSTRACT 

We report on the serendipitous discovery of a carbon star near the centre of the low- 
metallicity globular cluster NGC 6426. We determined its membership and chemical 
properties using medium-resolution spectra. The radial velocity of -159 km/s makes 
it a member of the cluster. We used photometric data from the literature and the 
COMARCS stellar atmospheric models to derive its luminosity, effective temperature, 
surface gravity, metallicity, and approximate C, N, and O abundance ratios. According 
to these properties, we suggest that this star is a genuine carbon rich low-metallicity 
AGB star. 
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ft Cen, RG055 and GRO70 jHardindll962l: iDickensll 1972ft . 
and a probable one in NGC 6402 (|C6te et all I1997T ). The 



two C H stars in ft Cen cannot be in the horizontal-branch 
phase (|Abiaet al.ll2003ft . Both GCs are of low-metallicity, 
and presumably belong to the Galactic halo. Sear ches for 
CH stars in other clusters have been unsuccessful (|Palmerl 
Il98d : iPalmer fc Whi3ll982ft . 

We report on the serendipitous discovery of a CH star 



in a globular cluster, and analyse its spectra to establish its 
properties and find clues as to its exact spectral type. 



2 OBSERVATIONS AND DATA REDUCTION 

The CH star is located about 24.5" or 1.6r c south-east of the 
centre of NGC 6426, a low-metallicity GC in the Galactic 
halo. Its coordinates, obtained with HST images and in the 
2MASS all-sky point source catalogu cQ, are: RA(J2000) = 
17' l 44 m 55 i :05; Dec(J2000)= +03°10!09'.'8. 

The star was discovered during an observing run on 10 
June 2010, at the 1.93m telescope of Observatoire de Haute- 
Provence (OHP) aimed at obtaining integrated spectra of 
Galactic globular clus ters with the CARELEC spectrograph 
|Lemaitre et al.lll990l ) and the grating 300 lines/mm. Two 
long-slit (5.5' x 2") spectra were obtained at the same posi- 
tion and PA=145.7°, both with 20-minute exposures and a 
seeing ~ 2.5". The corresponding dispersion and the spectral 
resolution were ~ 1.78A/pixels, and ~ 5A, and the spectral 
range was ~ 3700 — 6800A. Helium and Neon lamps were 
exposed at the beginning and the end of the night for wave- 
length calibration. For relative flux calibration and radial ve- 
locity calibration one flux standard HR8634 (|Hamuv et al.l 
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Figure 1. Spectrum of the CH star. The main spectral features 
are indicated. The spectrum of a classic CH star HD5223 from 
Goswami (2005) is shown in gray for comparison. 



I1992U1994T ), two Lick standard s t ars (H R5933 and HR7030) 
from the list of IWorthev et al] (|l994T), and an N-type CS 
(APM 2229+1902) from the list of iTotten fc Irwinl (|l99ct ) 
were observed on the same night. The signal-to-noise per 
resolution element at 5000A is S/N~29. 

In order to reach a better S/N in the blue and to benefit 
from a higher spectral resolution, additional long-slit spectra 
were obt ained with the Southern African Large Telesc ope 
(SALT) (|Bucklev et aLll200d ; lo 7 bonoghue D. et alj|2006r ) on 
12 and 2 April 2012 We used the Robert Stobie Spec- 
trograph (|Burgh et al] 120031 ') to make two 600 s. exposures 
with the grating GR900 and two 900 s. ones with GR2300. 
The slit width was 1.25" and the position angle PA=146°. 
The final reciprocal dispersions are 1.78 and 0.34 A pixel -1 
and the corresponding spectral resolution FWHM= 3 and 
1 A for the spectra taken with the two gratings, respec- 
tively. The reduction of the SALT long-slit data was done in 
the way described in lKniazev et al] |2008r ). Primary reduc- 
tion of the data was done with the SALT science pipeline 
l|Crawford et alj|201Ch . 

The data reduction and analysis of the OHP observa- 
tions were performed using MIDAS and IRAF. The disper- 
sion solution provides an accuracy of ~0.08 A for the wave- 
length calibration. Possible systematics in the wavelength 
calibration due to instrumental flexure were studied using 
the [OI] A5577 night sky line in the dispersion-corrected spec- 
tra, and by comparison of synthetic and observed spectra 
using the full spectrum fitting methods and constructing a 
det ailed line-spread fun ction of the spectrograph (LSF) (see 
e.g. iKoleva et aT]|2009r ). 

Finally the one-dimensional lower-resolution SALT 
spectrum was flux calibrated using the corresponding OHP 
spectrum and summed with it. The S/N in the resulting 
spectrum reaches ~15 at 4300 A, ~50 at 5000 A and ~160 
at 6300 A. The resolution of the spectrum obtained with the 
grating GR2300 is FWHM= lA and the mean S/N at 4740 
A is ~ 40. The final medium-resolution spectrum of the CS 
is shown in Fig[T] The high-resolution spectrum was used 
separately to estimate C 12 /C 13 , as explained in Sect. 5. 
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Figure 2. Radial velocity determination for the CH star using 
the method pPXF. 



3 RADIAL VELOCITY 

The heliocentric radial velocity of our star was derived 
by cross-correlation with the spectrum of the star APM 
2229+1902 using the iraf FXCOR package. Before cross- 
correlation we subtracted the continuum using the iraf 
task continuum and high-order polynomials. The cross- 
correlation shift (object-template) is 1.32 pixels. The radial 
velocity of our star is Vh = —163 ± 26 km/s, taking into ac- 
count a relative systematic shift (object-template) derived 
from fitting of the night sky lines (-1.2 - (-0.7)=-0.5 pix), and 
heliocentric corrections of -26 and -10 km/s, respectively, for 
our star and the template APM 2229+1902 (which has a ra- 
dial velocity Vh =-348 km/s). We also used a more robust 
pPXF method dCappellari fc Emselleml 120041 ) to measure 
the radial velocity on the same spectra. The result of the 
fitting is illustrated in Fig[2] Maximum penalized likelihood 
suppresses the effect of noise in the solution, so the accuracy 
is higher. Applying the aforementioned corrections, the re- 
sulting heliocentric radial velocity is Vh = —159 ± 5 km/s. 
Both estimates are close to t he radial velocity of NGC 6426 
V h = -162 + 23 l|Harrislll996l . 2010 edition). The probability 
that the CH star belongs to the field is very low, because 
using the Besangon model (|Robin et al.1120031 ). we find that 
the distribution of radial velocities of the 428 stars in the V 
magnitude range 10 to 18 located within 10' of the centre 
of the cluster is centered around a mean velocity of +13.22 
km/s, with a dispersion of 53.52 km/s. 



4 PHOTOMETRIC INFORMATION 

The photometric data for our CH star and the ones in NGC 
6402 and £1 Cen, a nd two classical Galact ic CH stars from 
the atlas of CSs bv lBarnbaum et al] (|l996l ) are listed in Ta- 
ble 1. The suc cessive lines are: (1) apparent visual distance 
modulus from lHarrisI (1 19961 . 2010 edition) (2) absolute mag- 
nitude in the V-band co rrected for Galactic e xtinction, (3) 
Galactic extinction from Schl egel et al] (|l998l ). (4-9) broad- 
band optical and infrared colors corrected for Galactic ex- 
tinction, (10) effective temperature. 

To estimate the coordinates of the CH star in NGC 
6402, we first obtained approximate coordinates on images 
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Table 1. Optical and 2MASS photometric data for CH stars. 
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N6426 


N6402 


RGO 55 


RGO 70 


HD5223 


V Ari 


(M - m) 


16.58 


16.71 


13.97 


13.97 


9.75 


9.19 


M v 


-2.58 


-1.86 


-2.39 


-2.36 


-1.43 


-1.14 


E(B-V) 


0.36 


0.60 


0.12 


0.12 


0.04 


0.14 


B-V 


1.69 


1.90 


1.56 


1.68 


1.43 


2.14 


V-I 


1.59 




1.45 


1.31 






V-K 


3.52 








2.82 


3.73 


J-H 


0.69 




0.67 


0.73 


0.53 


0.48 


H-K 


0.18 




0.23 


0.20 


0.17 


0.26 


J-K 


0.87 




0.90 


0.93 


0.70 


0.66 


Teff, K 


4100 




4468 


4287 


4360 


3600 



20 



from the CFHT archivtQ, and then identified it more accu- 
rately on HST images. These coordinates are: RA(J2000) = 
17 h 37 m 36 s 94 . Dec(J2000)= -03°14.'5'.'31. The colors and 
magnitudes f or stars in NGC 6426 i n the optic al bands were 
publis hed by lHatzidimitriou et al.l (1 19991 ) and iDotter et al.l 
(2011), using ground-based and HST images, respectively. 
Unfortunately, our star is satur ated on the HST images. In 
the color-magnitude diagram of lHatzidimitriou et all jl999l ) 
it is the reddest star in the cluster, slightly brighter than the 
stars at the tip of the red giant branch. Photometric data 
in the B, V, and I ban ds for the CH stars in NGC 6402 and 
Q Cen were taken from Cannon fe Stobid (|l973h .[ Cote et al.l 
(jl997h and IPancind (|2007h . Unfortunately, the CH star in 
NGC 6402 is too close to another star of same magnitude to 
obtain reliable near-infrared magnitudes from either 2MASS 
or Denial. The data for HD5223 and V Ari were extracted 
from the SIMBAD astronomical database. 

We calculated the near-infrared colors of the stars in 
the SAAO and TCS photometric systems using t he 2MASS 
colors and the transformation e q uation s from ICarpenterl 
|200ll ') and iRamfrez fc" Melendezl (|2004T ). With the colors 
thus derived, corrected for Galactic extinction, our star def- 
initely falls in the box delimiting CH stars in the two-color 
diagnostic diagram in the SAAO photometric system by 
iTotten et al.l (|200Cf) . Its effective t empe rature calculated us- 
ing relations from lAlonso et al.l fl999) and the (J — K) 
color in the TCS photometric system is T e tf= 4100 ± 12 5 
K. This corresponds to the class CH4 of iKeenanl {l993), 



which is equivalent to spectral type K4 III for oxygen stars. 
Similarly, the T ef f of RGO 55 and RGO 70 are 4468K and 
4287K, respectively. 



5 SPECTRAL ANALYSIS 
5.1 Classification 

The spectrum of the CH star was ex amined visually in terms 
of different spectral characteristics (|Goswamill2005f ) to avoid 
possible misclassification with C-N and R type CSs, which 
have spectra similar to those of CH stars. We compared 
the spectrum of the CH star with spectra at almost the 
same resolution for other such object s in the literature (e.g 
iBarnbaum et al . 1996; Goswami 20051) and found an approx- 
imate similarity in the shape of the main spectral features 
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Te=4000K, log(g)=0.5, C/0=10, N/0=0.5 
Te=4000K. log(g)=0.5, C/0=10, N/0=0.0 
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Figure 3. Observed medium-resolution spectrum of the CH star 
in NGC 6426 (black) fitted with the model one (green): T e ff = 
4000K, log(g) = 0.5, C/0= 10, N/0= +0.5. A model spectrum 
with lower N abundance is shown in red. 



with the spectra of two classical CH stars: HD5223 (type 
C-H3, C 2 index 4.5) and HD13826 (V Ari) (type C-H3.5, 
C 2 index 5.5). HD 5223 is a CH giant (see Fig. [[]) with 
[Fe/H]= -2.0 dex, C/0= 3.0, V Ari is a semi-regular pul- 
sating star with [F e/H] = -2.4, log(q) = -0.2, C/Q= 2, and 
12 C/ 13 C~ 6- 10 l|Aoki fc Tsuiilll997l : lGoswamill2005l ). The 
strength of the G band of CH (^4300 A) in the spectrum 
of our star, a mam characteristic feature of CH-type CSs, 
resembles that in the spectrum of HD 5223, and the same 
applies to the second branch near 4342 A. The line at 4226 
A is very weak. It is blended by molecular bands. The lines 
of atomic hydrogen and Ball are seen distinctly, which is not 
the case in C-R stars. The intensity of the CN band in our 
star is larger than in the case of HD 5223 and V Ari. The C2 
molecular bands near 4737, 5165, 5635, 6052 A are deeper 
than in HD 5223. This indicates a lower temperature for 
our star and /o r a higher C/O ratio l|Barnbaum et al1ll996l : 
lGoswamj|2005h . 

CH stars are classified into two types which follow dis- 
tinct evolutiona ry paths based on their 12 C/ 13 C ratios (e.g 
lGoswami|[2005h . Late-type CH stars with 12 C/ 13 C Js 100 
are intrinsic asymptotic branch (AGB) stars that produce 
the s-process elements internally. Early-type CH stars with 
low values 12 C"/ 13 C ^ 10 obtain the s-process elements via 
binary mass transfer. This is why the 12 C/ 13 C ratio is an 
important probe of stellar evolution. The band heads for 
12 C 12 C, 12 C 13 C" and 13 C 13 C* (at 4737A, 4744A and 4752A, 
respectively) are resolved even at medium resolution. In the 
following we will estimate 12 C/ 13 C and other parameters as 
accurately as possible using model atmospheres. 



5.2 Hydrostatic dust-free models 

To estimate the effective temperature, surface gravity, C/O 
and 12 C7 13 C for the CH star in NGC 6426 we used 
the hydrostat i c dus t-free models for carbon-rich giants of 
lAringer et al.l Ho09;) extended towards higher T e ff, log(g), 
and C/O to fit the observational data. We chose the fol- 
lowing model parameters: stellar mass IMq, metallicity 
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Figure 4. SALT spectrum of the CH star in NGC 6426 taken at 
the resolution of FWHM~ lA. The model with T eff = 4000K, 
l°9{9) = 0-5, C/0= 10, N/0= +0.5 is shown for comparison 
in grey. The band heads of 12 C 12 C, 12 C 13 C and 13 C 13 C are 
indicated by arrows. 



[Fe/H]= —2.0 dex, and solar abundance for the other 
elements includin g the isotopic ratio 12 C/ 13 C" = 86.8 
l|Scott et alJl2006h . The metallicity was chosen eq ual to the 
metallicity of NGC 6426 (given by lHarrislll996l . 2010 edi- 
tion). The model stellar mass roughly corresponds to the 
mass of stars at the Main-Sequen ce turnover point of NGC 
6426 (|Hatzidimitriou et al.l [l999h . Unfortunately, it is not 
possible to derive more accurate values of the mass and 
chemical abundance of the star, given the resolution of our 
spectra. 

The synthetic spectra computed with th e COMA code 
based on COM ARCS model atmospheres fsee lAringer et all 
I2009T ) were degraded by convolution with a Gaussian func- 
tion to fit the resolution of our data. The resulting synthetic 
colors were compared with the observed ones, providing ad- 
ditional arguments for a proper evaluation of the stellar pa- 
rameters. 

To obtain a good fit of the observational data with 
the model (both photometric and spectroscopic) we had to 
change not only the abundance of Carbon, but also those 
of N and O, because only increasing C at a given T B ff, 
log(g) did not allow us to fit the CN bands. A moderate 
increase of O with dre dge up or mass tran sfer is expected 
for carbon stars (e.g iLugaro et al.l |2012| . and references 
therein). Finally, we derived the following model parame- 
ters: T eff = 4000K, log(g) = 0.5, C/0= 10, [N/Fe]= +0.5 
dex, [0/Fe]= +0.5 dex which corresponds to an absolute 
magnitude My = —2.57 and the following synthetic colors: 
B — V = 1.71, V — K = 3.44, V - I = 1.57, J-H = 0.67, 
H-K = 0.24, J-K = 0.91. One can see from Table 1 that 
the agreement between the observed and model spectra and 
between the spectroscopic and photometric T e // has been 
reached in a wide spectral range. 

Fig. |3 shows how this model fits our medium-resolution 
spectra. Although our aim was not to estimate the abun- 
dances of different chemical elements on these spectra, the 
depth and shapes of the main molecular bands are fairly 
well adjusted. We did not vary the isotopic ratio 12 C/ 13 C, 
since there is no significant difference between the shape and 
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Figure 5. Survival of a binary star. The total encounter rate T 
is plotted vs the encounter rate 7 for single binaries. Globular 
clusters with concentration parameters c < 1.60 are indicated 
by triangles, the others by squares. The three globular clusters 
which harbour CH stars are indicated by full symbols. The units 
are arbitrary on both axes. 



depth of the C2 bands in the model and those in the spec- 
trum in the range 4737 - 4752 A. Fig. H shows the compari- 
son between the higher resolution spectrum of the CH star 
and the best fit model. The increased N abundance allows 
the strong CN features to be better matched. The derived 
luminosity, T e ff, log(g) and the high isotopic ratio 12 C/ 13 C 
indicate that this is probably a genuine AGB star. The de- 
rived T e f f i s normal for C-rich TP- AGB stars at a very low 
metallicity (|Marigo et al.ll200St ). 



6 SURVIVAL OF BINARIES IN GLOBULAR 
CLUSTERS 

We now examine the alternative possibility that the CH star 
is a binary star rather than a genuine AGB star. 

The dense environment of GCs favors the formation of 
binary stars, which in turn play an important role in their 
dyn amical evolution, by delaying the onset of core collapse 
fe.g. iHeggie et al1l2006t ). This dense environment may also 
lead to the dissolution of the wider binaries, so that the ac- 
tual number of bina ries depends on the central stellar den- 
sity and core radius (|VerbuntJ l2003), as well as on their sep- 
aration. 

In Fig. [5] we show the total encounter rate V oc p*' 5 ar 2 
versus the enc ounter rate 7 oc p°' 5 a/r c for single binaries 
(Verbunt 2003) for globular clusters in the halo (defined as 
having [Fe/H] < —0.80 dex). The parameters were calcu- 
lated u sing the core r adii r c and central luminosity densities 
p c from lHarrisI (Il996l, 2010 edition) and the mass-to-light ra- 
tios from lMcLaughlin fc van der Marell (|2005l ). We assumed 
that the M/L of NGC 6426, not available in the catalogs, 
is ~1.9, which is the most probable value for Galactic GCs 
at low metallicities. The uncertainty on this M/L is small 
(of the order of 10%). We also assumed that all binaries 
have the same semi-major axis a, which is a more drastic 
simplification. 

The clusters are aligned along a relation of almost unit 
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slope between the two encounter rates. The three GCs with 
known CH stars fall roughly in the middle of the relation. If 
all four CH stars are binaries, this suggests that no binaries 
can form when the total encounter rate is low, and that they 
are more rapidly disrupted than replenished if the encounter 
rate for single binaries is high. The number of known CH 
stars in clusters is admittedly too small to allow for any 
statistical inference on the probability of our CH star to be 
a binary. The three GCs have similar 7, but the parameter 
r is lower in NGC 6426 than in the two other clusters, thus 
making our star p erhaps les s likely to be a binary star. But, 
as pointed out by IVerbu nt (2003), r depends on the mass 
segregation, on the fraction of binaries and on their period 
distribution. Furthermore, both Y and 7 depend linearly on 
the semi-major axis, for which we assumed a universal value, 
but which can vary between a few and at least 200 Rq, so 
that presently a better way to determine the binary status 
of our star would be to monitor its radial velocity. 



7 CONCLUSION 

We have discovered a Carbon star of CH type in NGC 
6426 and derived its main physical and chemical parame- 
ters: M v = -2.58 mag, T efs = 4000-4100 K, log(g) ~ 0.5, 
[Fe/H]~ -2 dex, C/0~ 10, N/O- +0.5 dex, 12 C/ 13 C ~ 87. 
The data and the estimated encounter rates in the par- 
ent GC indicate that the object is likely an intrinsic low- 
met allicity carbon-rich AGB star, but additional exten- 
sive high-resolution observational and theoretical studies are 
needed to reveal the role of a possible presently invisible 
companion on its evolution. 
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